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Abstract 


A method for the on-line measurement of the heat of combustion of gaseous 
hydrocarbon fuel mixtures has been developed and tested . The method involves com- 
bustion of a test gas with a measured quantity of air to achieve a preset concentration 
of oxygen in the combustion products . This method involves using a controller which 
maintains the fuel (gas) volumetric flow rate at a level consistent with the desired oxy- 
gen concentration in the combustion products. The heat of combustion is determined 
from a known correlation with the fuel flow rate. An on-line computer accesses the 
fuel flow data and displays the heat of combustion measurement at desired time inter- 
vals. This technique appears to be especially applicable for measuring heats of com- 
bustion of hydrocarbon mixtures of unknown composition such as natural gas. 


Introduction 

Several well-known techniques are available cur- 
rently for the measurement of heat of combustion of 
hydrocarbon fuels. Constant-volume bomb calorimeters 
are commonly used for determining heat of combus- 
tion of liquid fuels (refs. 1 and 2). Similarly, constant- 
pressure flame calorimeters are employed for measuring 
the heat of combustion of gaseous fuels (refs. 3 and 4). 
Singh et al. (ref. 5) have developed an on-line technique 
which requires injection of pure oxygen prior to combus- 
tion for measuring the heat of combustion of gaseous 
fuels. The injection is necessary to maintain the required 
21 -percent oxygen concentration in the combustion 
products. Singh et al. have also developed a simplified 
method that eliminates the need for the addition of 
pure oxygen and uses a computational algorithm that is 
valid only for a prescribed air flow rate and combustion- 
product oxygen concentration (ref. 6). The present tech- 
nique offers an improvement over this latter method by 
using a computational algorithm that is valid over a wide 
range of flow rates and oxygen concentrations. Develop- 
ment of the present measurement system satisfies the 
need for accurate, on-line measurement of heat of com- 
bustion of gaseous fuel mixtures with methane as the 
dominant component. This fuel has been used in the 
combustor of the Langley 8 Foot High-Temperature 
Tunnel (8-Ft HTT) to produce temperature levels 
approaching 2000 K. 

In the present study, an on-line gas calorimeter was 
developed and tested for high-precision measurement of 
heat of combustion of gaseous hydrocarbon mixtures. 
The method described herein involves the measurement 
of the concentration of oxygen in the combustion prod- 
ucts and its correlation with the heat of combustion of the 
fuel. 

Symbols 

m m 

A = 4£M + (l+P)2M-4i > 

i - 1 i= 1 


Cf fuel flow rate correction factor 

CyH y hydrocarbon 

E optical encoder data 

/ total hydrocarbon mole fraction 

fl mole fraction of ith hydrocarbon species 

/ inert fraction of fuel 

/ air volumetric flow rate 

m number of hydrocarbon species 

n f fuel volumetric flow rate 

P combustion-product oxygen mole fraction 

P 0 mole fraction of oxygen in air 

R 2 coefficient of determination 

a, P regression coefficients 

AH heat of combustion 

Experimental Setup 

Die prototype gas calorimeter is shown schemati- 
cally in figure 1 and photographically in figure 2. The 
fuel gas is mixed with air and fed into a specially 
designed burner. The fuel flow rate is measured by a vol- 
umetric flowmeter with an optical encoder mounted on 
its shaft. The flow rate is determined by measuring the 
counts/min of the etched lines on the encoder disk. Each 
revolution of the shaft, which is nominally equivalent to 
1000 cc, produces 1200 counts. 

A burner that is fabricated with an internal honey- 
comb structure was used to ensure a stable flame over a 
range of air and fuel flow rates. The burner was enclosed 
in a vented glass container which prevented atmospheric 
contamination of the combustion products. 

The oxygen concentration in the combustion prod- 
ucts was measured by pumping the combustion products 
through a yttrium-stabilized zirconium oxide oxygen 
sensor that produces a voltage related to the concentra- 
tion of oxygen in the environment to which it is exposed 
(ref. 7). The oxygen sensor output was directed to a 



controller that maintained the desired concentration of 
oxygen in the combustion products by varying the fuel 
flow rate. This technique requires a sufficient quantity of 
oxygen (>0.1 percent) in the combustion products to be 
measurable by the oxygen sensor. A set point concentra- 
tion of 3-percent oxygen was used throughout this study. 

Because of the precision uncertainties of the volu- 
metric flowmeter, the oxygen sensor, and the air flow- 
meter, and because of the inherent unsteadiness of the 
combustion process, the oxygen concentration in the 
product gases exhibited small fluctuations about the 
3-percent set point. To account for the random fluctua- 
tions in the oxygen concentration, the mean and the stan- 
dard deviation (of heat of combustion) were obtained for 
a set of data points at several time intervals. 

Experimental Procedure 

The oxygen sensor was calibrated with mixtures of 
known oxygen concentrations so that a sensor output 
(in mV) versus oxygen-concentration calibration curve 
could be obtained. Subsequently, the system was flushed 
with air for several minutes to scavenge any fuel that was 
left inside the system from the previous experiment. Air 
from a pressurized cylinder was used for combustion, 
and its oxygen concentration was measured with the oxy- 
gen sensor prior to the experiment. The fuel flowmeter 
was also calibrated so that the number of counts/min 
from the optical encoder could be converted readily into 
flow rate in standard cubic centimeters per minute 
(seem). 

A given air flow rate from the pressurized cylinder 
was initiated, and the fuel flow rate was gradually 
increased to a level at which a stable flame was estab- 
lished. The 3-percent oxygen concentration was dialed as 
a set point in the controller, and the controller was acti- 
vated to maintain the fuel flow rate. The fuel flow-rate 
data in counts/min were recorded at regular intervals 
and were used in a correlation (described in the “Test 
Results” section) to obtain the value of heat of combus- 
tion. Over a time period, several data points were 
recorded, and the mean and the standard deviation were 
computed for the period. 


where / is the total hydrocarbon mole fraction, m is the 
number of hydrocarbon species,/* is the mole fraction of 
the ith hydrocarbon species (£/* =1), and A H t is the heat 
of combustion of the ith hydrocarbon species. 

The combustion of this mixture with air 1 may be 
expressed as 

/ + n[f{ fl C x H yi + f 2 C x H y2 + ... •+ f m C x H y ) + (1 -/)/] 

Heat 

«/(/ 1*1 + f 2 x 2 f m x m )CO 2 

+ f(/ I y,+/2J-2+-+/ m J’JH 2 o 

+ [ lP o~ »/(/ 1*1 + • • ■ + fm x m) - jif iyi + - + 

+ [/(i-p 0 )iN 2 +n(i-/)/ (2) 

where / is the volumetric flow rate of air, n is the volu- 
metric flow rate of the gaseous fuel mixture, / is the non- 
combustible inert fraction of the test gas, and P 0 is the 
mole fraction of oxygen in the air. 

The mole fraction of oxygen in the combustion prod- 
ucts P can be expressed as 


4lP o - An f X f i x i -nf'lL f i y i 


P = 


1 = 1 


1 = 1 


«/XM + 4»-4n/ + 4! 
i=l 


Solving for/, 
/ = 


4l(P g -P)-4nP 


4 

L j=l f=l 


(3) 


(4) 


Substitution for / in equation (1) yields the following 
expression: 








4 X f i X ‘ + ( 1 + X f i y i ~ AP 


l ;=i 


(5) 


Theoretical Principles 

The heat of combustion A H^x a gaseous mixture 

that has several saturated hydrocarbons and a noncom- 
bustible inert gas may be expressed as 

m 

AW mix = f X f i AH i (D 

/=! 


Let 

m m 

A = 4^ ffy + (1 +P) £ f i y i - 4 P 
i=l i=l 

A is a measure of the effective hydrocarbon content of 
the combustible component. Note that for saturated 


1 Air is assumed to be composed of nitrogen and oxygen. 
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hydrocarbons, A can be further simplified since 
y t = 2x, + 2. 


Table I. Values of AH, A, and AH I A for Pure 
Saturated Hydrocarbons 


The Affjnjx may now be expressed as 

A1,n ™ Y f i AH i 
2, A~ 

i=l 


mix = 


r^ K r 0 -rj 


-4 P 


( 6 ) 


For a binary mixture of a single pure hydrocarbon 
species (m = 1) and an inert gas, the above equation 
reduces to 


AH 


mix 


r 4 KPp-P) 

n 



(7) 


where A = 4x ■+ (1 + P)y - 4P, and AH is the heat of com- 
bustion of the hydrocarbon species. Note that for meth- 
ane (y = 4), A is independent of P. 

For a mixture of unknown composition, the value of 
X, ffAHi /A to be used in equation (6) is not known since f t 
and A are unknown. Table I shows values of AH (ref. 8), 
A, and AH /A for various saturated hydrocarbons. (When 
more than one form of a hydrocarbon exists, such as 
n-butane and isobutane, the average value of AH is 
given.) A value of 0.03 for P is used in the calculation of 
A in each case. In the present technique, the adoption of a 
universal value of the X/*AH;/A coefficient in equa- 
tion (6) is crucial for measuring the heat of combustion 
of a gaseous mixture of unknown composition. However, 
as noted from table!, AH/A varies from 25.94 for hexane 
to 26.6 for methane. For the present technique it is suffi- 
cient to adopt a value of this coefficient that will mini- 
mize the error in the prediction of the A H^ for the 
domain of samples to be tested. Thus, for the domain of 
samples where all species occur with equal probability, 
the mean value of 26.27 is the appropriate choice for 
' 'LfiAHj/A in equation (6). 

In nature, however, hydrocarbons are not distributed 
uniformly. Methane makes up the bulk of the hydro- 
carbon component in most natural gas samples. Figure 3 


Hydrocarbon 

AH* 

kcal/mol 

A 

(P = 0.03) 

AH/A 

Methane, CH 4 

212.80 

8 

26.60 

Ethane, C 2 H 6 

372.82 

14.06 

26.51 

Propane, C 3 H 8 

530.61 

20.12 

26.37 

Butane, C 4 H 10 

686.65 

26.18 

26.23 

Pentane, C 5 H 12 

837.05 

32.24 

25.96 

Hexane, C 6 H 14 

993.51 

38.30 

25.94 




Mean = 26.27 ±0.51* 


*A11 heats of combustion AH are given as the gross heat of combustion 
(producing liquid water at 25°C (ref. 8)). 

*A11 errors in this paper are given at the 95-percent level of confidence. 


shows the variation of X/-AH; with the methane mole 
fraction of the hydrocarbon component for 265 well sam- 
ples that were analyzed in a 1982 Bureau of Mines sur- 
vey of natural gases (ref. 9). Figure 4 shows the variation 
of ZffiHi/A for these samples. Although the value of 
X fjAHi varies significantly from sample to sample, the 
value of XyjAH; /A varies in a narrower range. Figure 5 
presents a frequency distribution plot of Xj^A H t /A for 
these samples. More than 90 percent of these well 
samples have a 'Zf i AH i /A between 26.5 and 26.6. Thus, 
the mean value of 26.55 is the appropriate choice for 
X, ffAHi/A for a wide range of natural gas samples. 

Table II shows values of X^AH), A , and X^-AH; /A 
for six natural gas pipeline samples from the Bureau of 
Mines survey. Although values of X^AH^andA vary sig- 
nificantly from sample to sample, the value of X/*AH* /A 
is constricted to a narrow band from 26.53 to 26.59, 
in agreement with the mean value cited previously. 
Thus, although the adoption of a universal value for the 
X/jAH/A coefficient inevitably introduces an error in the 
calculated value of AH, this error will be negligible in the 
overwhelming majority of samples from natural gas 
sources. 


Table 31 The 'LfjAHj, A, and X/Afy /A for Six Natural Gas Pipeline Samples 


Hydrocarbon component concentrations,* mole percent 

I/iAffj, 

kcal/mol 

A 

(P = 0.03) 

XfiAHj/A 

No. 

ch 4 

c 2 h 6 

c 3 h 8 

c 4 h 10 

C 5 H 12 

c 6 h 14 

1 

95! 

3! 

0.5 

0.2 

0 ! 

0.1 

219.9 

8.35 

26.59 

2 

86.5 

8.2 

2.7 

1.0 



240.4 

9.15 

26.55 

3 

86.6 

8 ! 

2.6 

1.0 

0.5 

0.2 

241.3 

9.18 

26.55 

4 

81.8 

10.0 

4.1 

1.8 

0.8 

0.3 

255.2 

9.74 

26.53 

5 

89.6 

6.5 

1.8 

0.7 

0.3 

0.2 

233.8 

8.88 

26.56 

6 

93.6 

3.7 

0.6 

0.3 

0 ! 

Trace 

219.2 

8.39 

26.59 










Mean = 26.56 ± 0.04 


*From reference 9. 
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Test Results 

In the present study, l = 7000 seem and P = 0.03, For 
these values equation (6) reduces to 


A H 


mix 


r 28000(P - 0.03) 1 jn f AH; 

[ i — <8) 

1 = 1 


Equation (8) forms the basis for calculation of the A H mix 
from the volumetric fuel flow rate measurements. 

The present technique was first applied to two pure 
gas samples: methane and ethane. Since heats of com- 
bustion of these gases were already known, the results 
were used to determine the precision error of the method, 
and to this end a correction factor Cy was introduced in 
the calculation of volumetric flow rate from the optical 
encoder data E to correct for bias in the system. Table HI 
shows the experimental data for methane that involves 
optical encoder data E in counts/3 min intervals, the 
fuel flow rate n in seem calculated from the optical 
encoder data (n = E/3 min * 1000 cc/1200 * Cj), and the 
heat of combustion (A H) in keal/mol calculated from 
equation (8); P Q was measured at 0.209. The mean value 


Table HI. Methane Data 


E , counts/ 3 min 

Fuel flow rate, 
n* seem 

AH, kcal/mol 

2125 

621.4 

211.7 

2118 

619.3 

212.4 

2108 

616.4 

213.5 

2115 

618.4 

212.7 

2116 

618.7 

212.6 

2108 

616.4 

213.5 

2110 

617.0 

213.3 

2118 

619.3 

212.4 

2112 

617.6 

213.1 


Mean = 618.3 ±3.1 

Mean = 212.80 ±1.22 


*Cf= 1 . 053 . 


of the heat-of-combustion data is 212.80 ±1.22 for 
methane. 

Table IV shows the data for ethane; P 0 was measured 
at 0.210. Readings were taken at 5-min intervals. The 
mean value of the heat-of-combustion data for ethane is 
372.82 ± 2.60, These two cases indicate that a precision 
approaching 0.6 percent can be achieved for pure satu- 
rated hydrocarbon fuels such as methane and ethane. 


Table IV. Ethane Data 


E, counts/ 5 min 

Fuel flow rate, 
n * seem 

AH, kcal/mol 

2122 

356.2 

370.9 

2119 

355.7 

371.4 

2121 

356.0 

371.1 

2122 

356.2 

370.9 

2117 

355.4 

371.8 

2112 

354.5 

372.7 

2113 

354.7 

372.5 

2112 

354.5 

372.7 

2107 

353.7 

373.5 

2102 

352.8 

374.4 

2108 

353.8 

373.4 

2103 

353.0 

374.3 

2103 

353.0 

374.3 

2101 

352.7 

374.6 

2104 

353.2 

374.1 


Mean = 354.4 ±2.5 

Mean = 372.82 ±2.60 


*C f = 1 . 007 . 


The system was also used to find the heat of combus- 
tion of a mixture that was 92.49-percent methane and 
7.51-percent nitrogen, which has a calculated heat of 
combustion (eq. (1)) of 196.82 kcal/mol; P Q was mea- 
sured at 0.209. Table V shows data for this mixture, 
along with the fuel mole fraction / calculated from 


Table V. Data for 92.49 Percent CH 4 Plus 7.51 Percent N 2 


E , counts/ 3 min 

Fuel flow rate, 
n* seem 

AH, kcal/mol 

Fuel mole fraction,/ 
(calculated from eq. (4)) 

2297 

672.8 

195.0 

0.9162 

2265 

663.4 

197.8 

0.9294 

2284 

669.0 

196.1 

0.9215 

2288 

670.1 

195.8 

0.9199 

2260 

661.9 

198.2 

0.9315 

2246 

657.8 

199.5 

0.9374 

2262 

662.5 

198.0 

0.9306 

2301 

673.9 

194.6 

0.9146 

2280 

667.8 

196.4 

0.9231 


Mean = 666.6 ± 10.3 

Mean = 196.82 ± 3.10 

Mean = 0.9249 ±0.0145 


*C f = 1 . 054 . 
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Table VI. Summary of Data for Methane, Ethane, and Methane Plus Nitrogen 


Atf est = 48.6 + 


3.42 x 10 J 


Regression output: a, 48.5844; Standard error of Ai7 est , 2.005064; R 2 , 0.999463; 
Number of observations, 33; Degrees of freedom, 31; p, 342320.7; Standard 
error of P, 1425.077 


Calibration gas 

Flow data (P), 
counts/3 min 

AH, kcal/mol 

Atf est , kcal/mol 

Ethane 

1061.0 

372.82 

371.22 



1059.5 

372.82 

371.68 



1060.5 

372.82 

371.38 



1061.0 

372.82 

371.22 



1058.5 

372.82 

371.99 



1056.0 

372.82 

372.75 



1056.5 

372.82 

372.60 



1056.0 

372.82 

372.75 



1053.5 

372.82 

373.52 



1051.0 

372.82 

374.29 




372.82 

373.37 




372.82 

374.14 



1051.5 

372.82 

374.14 

\ 

/ 

HHHDSSsSflMH 

372.82 

374.45 


/ 


372.82 

373.98 

Methane 

2125.0 

212.80 

209.68 




212.80 

210.21 



2108.0 

212.80 

210.98 



2115.0 

212.80 

210.44 



2116.0 

212.80 

210.36 



2108.0 

212.80 

210.98 



2110.0 

212.80 

210.82 

\ 

/ 

2118.0 

212.80 




2112.0 

212.80 


1 Methane, 92.49 percent 

2297.0 

196.82 

197.61 



2265.0 

196.82 

199.72 



2284.0 

196.82 

198.46 



2288.0 

196.82 





196.82 




2246.0 

196.82 

201.00 



2262.0 

196.82 

199.92 

\ 

/ 

2301.0 

196.82 

197.35 

2280.0 

196.82 

198.73 


equation (4). A value of 212.80 for AH and a value of 8 
for A were used since methane was the only hydrocarbon 
involved. The mean value of the heat-of-combustion data 
is 196.82 ± 3.10. The mean value of the fuel mole frac- 
tion data is 0.9249 ± 0.0145. 

Numerical Approach 

The heat of combustion may also be expressed as 

AH = a + 1 (9) 

E 


where E is flow data, either in the form of flow units or 
directly in optical encoder counts, and a and P are 
constants determined from a least-squares regression 
analysis of flow data and heat of combustion for gases of 
various heat content, assuming constant values for /, P 0 , 
and P. Equation (9) allows calibration of AH directly in 
terms of optical encoder counts using calorimetric gas 
standards, thus obviating the intermediate steps of quan- 
tifying the flow rate and oxygen concentration of the 
combustion air, calibrating fuel flow rate from optical 
encoder counts, and calibrating oxygen concentration 
from oxygen sensor voltage. This regression algorithm 
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was applied to optical encoder readings for the methane, 
ethane, and methane-plus-nitrogen data presented above 
and effectively provides a 3-point calibration of AH in 
terms of E (optical encoder data for ethane were normal- 
ized for 3-min intervals). Table VI summarizes these 
data and presents results of the regression analysis. A 
graphical plot of the data, the least-squares regression 
curve fit AH est , and the ±2 standard error of A H cst (SEE) 
band, which approximates the 95-percent confidence 
interval about the curve fit (ref. 10), are shown in 
figure 6. 

Conclusions 

An experimental technique has been developed and 
tested to measure the heat of combustion of a saturated 
hydrocarbon and inert gas mixture. The technique is 
especially attractive for applications in which an on-line 
measurement of the heat of combustion of natural gas is 
required. The system worked reasonably well during the 
course of laboratory investigation and thus validated the 
technique. Preliminary test results indicate that this tech- 
nique can determine the heat of combustion of a gas of 
unknown composition with an uncertainty on the order of 
2 percent. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
January 22, 1996 
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Figure 1. System diagram. 
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Figure 2. Experimental setup. 
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Figure 3. Variation in calculated heat of combustion IfiAHj with methane mole fraction of hydrocarbon component for 
265 natural gas well samples. 
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Figure 5. Distribution of IfjAHj/A for 265 natural gas well samples. Mean = 26.55; Standard deviation = 0.06. 
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